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Abstract 

Statistical properties of turbulent magnetic fields in radio-synchrotron sources should imprint on the statistics of polarimetric ob- 
t-H servables. In search of these imprints, we calculate correlation and cross-correlation functions from a set of observables containing 

t-H the total intensity /, the polarized intensity P and the Faraday depth tp. The correlation functions are evaluated for all combinations 

of observables up to fourth order in the magnetic field B. We derive these as far as possible analytically and from first principles 
only using some basic assumptions such as Gaussian statistics of the underlying magnetic field in the observed region and statistical 
homogeneity. We further assume some simplifications to reduce the complexity of the calculations, as for a start we were interested in 
a proof of concept. Using this statistical approach, we show that it is in principle possible to gain information about the helical part of 
the magnetic power spectrum, namely via the correlation functions (P(]i ± )<p(k' ± )<p(k"))v and (/(k ± )0(k^)0(k")) B . Using this insight, 
we construct an easy-to-use test for helicity, called LITMUS (Local Inference Test for Magnetic fields which Uncovers heliceS) which 
gives a spectrally integrated measure of helicity. For now, all calculations are given in a Faraday-free case, but set up in a way so that 
Faraday rotational effects could be included later on. 

Key words. Magnetic fields - Methods: data analysis - Methods: statistical 

o | 1. Introduction 

q Magnetic fields are observed in almost all astronomical objects, they permeate planets and stars as well as galaxies and clusters. 
5— i Most, if not all, of the interstellar and intergalactic plasma appears to be magnetized and the magnetic fields contribute significantly 
to physical processes. Examples include the formation of stars ( Price et al.|2009| l, the anisotropy of transport processes (thermal 



d conduction or plasma resistivity, see e.g. Narayan & Medvedev| ( 2001| )), the angular momentum transport in accretion discs or the 
propagation of cosmic ray populations ( Strong et al. 2007][ 

Although magnetic fields are ubiquitous in the cosmos, we often cannot treat them properly in astrophysical situations due to 
h> the lack of knowledge of their properties. Cosmic magnetic fields are difficult to observe and their distribution, evolution and origins 
are far from being perfectly understood. We have three main sources of information: the Zeeman effect, synchrotron radiation and 
Faraday rotation. The Zeeman effect is extremely difficult to detect, because other line shifting effects, such as thermal Doppler- 
C^l broadening, are usually stronger. We obtain a great deal of information from synchrotron radiation but only regarding the magnetic 
field component perpendicular to the line of sight. In order to get a picture of the 3D magnetic field, one needs another source of 
qq information. This leads us to Faraday rotation, the change of the polarisation-plane of long wavelength radiation due to a magnetic 
field along the line of sight. Faraday rotation provides a powerful tool, but is also difficult to observe, to evaluate, and to interpret 
f^") due to the involved line of sight projection. This projection is one of the main obstacles to understand the 3D properties of magnetic 
i-H fields. 

One important property of cosmic magnetic fields that we do not know much about is magnetic helicity. It is defined as the 
. integral 

X ^ r 

^ //= J A-Bdx 3 (1) 

over a Volume V with surface dV on which n ■ B = 0; where A refers to the vector potential from electrodynamics with B = V X A. 
Helicity is a measure for the "spiral quality" of a magnetic field. It quantifies how much the magnetic field lines are sheared and 
twisted and counts the number of spirals the field lines exhibit within a given volume. Particularly turbulent magnetic fields should 
show considerable helicity. The relevance of helicity has increased since its inclusion as an essential element in the magnetic dynamo 
theory, which tries to explain the sustainement of magnetic fields on large scales over cosmic timescales (see Subrama nian|2002"t 
Brandenburg & Subramanian 2005a b). The possible operation of a large scale dynamo for instance is directly connected to the 



generation of helicity in turbulent environments (see Shukurov et al. 2006; Brandenburg 2009; Sokoloff 2007 ), a process which, 
until now, could not be verified through observation. 

This study is particularly concerned with the question of how to extract knowledge regarding turbulent magnetic helicity spectra 
from the statistical information found in radio-observational data involving polarisation and Faraday rotation measurements. It was 
highly motivated by the studies of Volegova & Stepanov (2010). 
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Information on magnetic fields can be imprinted onto radio data by two of the processes already mentioned above: Synchrotron 
emission and Faraday rotation. Yet information is not only contained in their mean values but also in higher order correlation and 
cross-correlation functions. 

We therefore investigate a set of suitable radio observables for their cross-correlations, to see how these are connected to the 
statistical properties of the magnetic fields to be examined. This idea goes back to previous works by Spangler] ( 1982 1983); 
Eliek| ( |1989a|b| l; |Enfilin & Vogt| p003"l >; |Kahniashvili & Vachaspatrj ( |2"00r5] i; |Waelkens etal] ( |2009| ). The set of radio observables 



we investigate contains the total intensity I(x), the polarised intensity P(x) and the Faraday depth <p{x). We work out all correlation 
functions between them in a general framework. We restrict ourselves to fourth order in the magnetic field strength and as far as 
possible we do all calculations analytically. The aim is to find a direct relation to statistical properties of the magnetic fields, such 
as their power spectra. 

The intensity I(x) and the polarised intensity P{x) are connected to the synchrotron emission within a magnetized volume. We 
assume them to be taken at sufficiently high frequencies and, therefore, free of Faraday rotation. The Faraday depth <p(x) is measured 
via the Faraday rotation of a polarized background source at a different frequency seen through the same volume. The observational 
situation is visualized in Fig. Q, 

With regard to these observable quantities, we can successfully establish all correlation functions in the form of analytical 
relations to the magnetic field power and helicity spectra implementing Gaussian field statistics for simplicity. The result here is to 
prove that it is possible, in principle, to gain information not only in respect of the total but also regarding the helical part of the 
magnetic power spectrum, namely via (P(k ± )0(k^)</>(k"))B and (7(k ± )0(k^)0(k"))B- 

Gaussian magnetic fields statistics is not what numerical simulations of MHD turbulence find (see Waelkens et al. 2009). 



However, they are the starting point of any analysis of high-order correlation functions. In case all we know statistically about the 
fields is their two point correlation, the only assumption which expresses solely this knowledge is a Gaussian with the correlation 
tensor being the covariance matrix. Any other distribution function would contain more information in a Shannon-Boltzmann sense. 

If additional information on higher order statistics is available, this could be incorporated via perturbative methods. These would 
expand around the Gaussian case, which has therefore to be worked out first, as we do in this work. 

Based on our results we further present the LITMUS test (Local Inference Test for Magnetic fields which Uncovers heliceS), 
a first simple procedure to probe data for helicity. An analysis of real and simulated data using this test along with a thourogh 
investigation of its applicability can be found in Oppermann et al. (20101. The study is organised as follows: Section [2] presents 
our method and the general formalism we developed to evaluate the correlation functions analytically. Section|3]details a complete 

example calculation for one of the correlation functions, namely (P(k ± )-P(k^))£. SectionJ^Jthen presents all the correlation functions 
up to fourth order in magnetic field strength. Section|5]introduces the LITMUS test. Section|6]presents finally a thorough evaluation 
of our findings. Details of the derivation of the other correlation functions are listed in the Appendix which also contains the 
remaining technical information regarding the study. 




Figure 1. Schematic of an observational situation for which our set of correlation functions is suitable (modified picture taken from 
Wael kens et"aT1p009] l). 



2. Methods 

We now proceed to calculate the correlation functions of I(x ± ), P(x ± ) and <p(x ± ). Since all calculations resemble each other in 
respect of certain basic features, a general framework has been developed for them. Before presenting that, we introduce our basic 
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notations and the magnetic correlation tensor, a quantity that will be referred to frequently but requires some preliminary explanation 
due to its complexity. 



2.1. Notation 

During this study we use the following definitions for the radio observables 7(x ± ), P(x ± ) and <p(x ± ): 



dz {B\ + B\), (2) 
= j dz{B\-Bl+2iB x B 2 ), and (3) 



= dzB 3 . (4) 



Throughout this study, the coordinate axis z always equals the line of sight. For convenience, all fore factors are suppressed including 
the electron density n e , assumed to be constant. A detailed derivation can be found in Appendix [B| 

Apart from this, we like to introduce some further notation we use frequently. In the following, the vectors r or r' shall always 
denote a combination such as x' - x to be defined exactly when needed. Furthermore, w = (k"',0), u = (k",0), v = (k^,0) and 
a = (-q ± - k^, -q z ) holds. 



2.2. The magnetic correlation tensor 

We assume generally Gaussian statistics for the magnetic field distribution in the observed region. This is, of course, a simplification 
but represents an initial starting point, especially as we are mainly interested in a proof of concept and therefore seek for simplicity. 
In addition, we assume statistical homogeneity to first establish the most illuminating cases. This assumption is widely used in the 
literature. For an arbitrary field iff, statistical homogeneity means that the two point correlation function of the field depends only on 
the distance of the two parts, (tfr(x')il/*(x)) = C(r) with r = x' - x. This automatically implies for this correlation function in Fourier 
space: 

<<«k>*(k)> = (270 3 5 3 (k' - k)/y(A0 



where the (/(-power spectrum is specified by the Fourier transformed correlation function P^ik') oc J dr'C^(r) exp [ik'r] as stated by 
the Wiener-Kinchin-theorem (Bracewell 2000]l. 

Within this study, the magnetic correlation tensor M, ; (x, x') = (Bj(x)Bj(x')) is frequently used for which this translational 
invariance leads to 

Mtj(x,x') = Mij(x' — X = r) in normal space, and (5) 



My(k, k') = (27r) 3 5 3 (k' - k)My(k') in Fourier space. (6) 
For homogeneous and isotropic magnetic turbulence the translationally invariant magnetic correlation tensor can be written as 

Mij(r) = M N (r)6 u + (M L (r) - M N (r)) ^ + M H (r)e ijm r m (7) 

with the longitudinal, normal and helical spectra denoted by Mi(r), M^(r) and Mn{r) respectively. The solenoidal condition V B = 
enables the connection of the two non-helical spectra by M N (r) = j- ^(r 2 M L (r)). By applying a Fourier transformation, we obtain: 

My(k) = M N (k){5 u - ~y~) - ie ijm H(k) ^ . (8) 

In this case, the condition V ■ B = was used directly in the form kjMjj(k) = to reduce the degrees of freedom to the normal and 
the helical spectra. These two functions are specified in terms of their real space counterparts as 

M N (k) = Jdr 3 M N (r) exp[/kr], and (9) 
" (k) = dk^ H(k) =d~k\ driMll{r) ex P[ /kr ]- ( 10 ) 



Some interesting properties are: 



Mjv(O) = M L (0) by definition, (11) 
M y (0) = M N (0) Sjj, and (12) 
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Mn(0) = 0, since would diverge otherwise. (13) 

The magnetic correlation tensor is closely related to the energy spectrum of the magnetic field. The field's mean energy density can 
be expressed as follows 

s< B2 «> - c / <&< Z - c / H Z *«*> = 

= 7T^3 dkk 2 M N (k)= dke B (k). (14) 

(27TP Jo Jo 

In this case, we have used Mu(k) = 2M^ and e#(fc) denotes the ID-energy density of B. From this we derive 

k 2 M N (k) 

£b() = &n 3 ' ( } 

which is used in the following to replace M^ik) by the more commonly applied magnetic energy spectrum. 

Analogously to ( 14 1, we can relate the helical part of the spectrum H{k) to the mean current helicity j -B and deduce a ID-he lical 
energy density enik)' 



(j ■ B) = (B • (V X B)> = (BMeujdrfifr + r)>| r=0 = 3,, e/ , 7 M y (r)|, =0 
, uk 3 r ^ 3 

= «iy 



kk 



/dk f dk I 

p^kiMijQQ = J —^nje ljm H(ky 

I 

Jo 



IStm 

J 

dke H (k). (16) 



We can read off the ID-helical energy density e H {k) which also can be used to substitute H(k): 

e H (k) = j— . (17) 

Nevertheless, in our calculations we want to relate to the helicity spectrum R(k) rather then to the current helicity spectrum H(k) 
since it is the helicity B ■ A that usually is the subject of magnetohydrodynamics and dynamo theory. 

Fortunately, for isotropic, turbulent fields the current helicity and the magnetic helicity B ■ A are closely connected. Since the 
current helicity (V x B) • B has the same mathematical structure as the helicity (V x A) • A, we can perform exactly the same derivation 
as in ( 16 1 to show that the mean helicity relates to the helicity spectrum R(k) of the correlation tensor for the vector potential in the 



same way as the mean current helicity relates to the current helicity spectrum of the magnetic correlation tensor. We can construct 
the correlation tensor (A(x)A*(x')) for the magnetic vector potential A similar to the magnetic correlation tensor. If we assume 
translational invariance for the statistics and apply the Lorenz gauge condition V ■ A = to the vector potential, all deductions made 
for the magnetic correlation tensor will also hold for the correlation tensor of the vector potential. In Fourier space it will have the 
form: 

(A m (k)A;(k')> = R mn (k) = R N (k){5 mn - - ie mm ,R H (k)j (18) 

From here we can start by rewriting the magnetic correlation tensor in terms of the correlation tensor for the vector potential: 

<B,(x)B}(x')> = e ilm e jrn d Xl d Xr (A m (xW n (x')) (19) 

If we now perfom a Fourier transformation on both sides and furthermore apply the condition for translational invariance (|6]l we 
get: 

/dk 3 C dk 3 

— — My(k)= — RmniV. (20) 

(Zny J (Zny 

Thus, we find 

M/y(k) = e ilm £j m kik r R mn (k). (21) 

Using 



we can rewrite the right hand side of (21 



k k 

£u m £j m kik r R mn (k) = £ii m £j rn kik r (6 m „ - -p-) - ki£ jr ik r k i R H (k)/k + ik 2 £ ijr k r R H (k)/k (23) 

=0 antisymmetric in i,j 



symmetric m i,j 
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Finally, we find the relation between the current helicity spectrum H(k) and the helicity spectrum Rh by identifying the antisym- 
metric part of the left hand side of (21 1 with the antisymmetric part of the right hand side: 



-ie ijm k m H(k)/k = +ik ei jr k r R H (k)/k 
=> H(k) = -k 2 R H (k) 



(24) 



From now on, although we mostly speak of helicity, we actually deal with current helicity for convenience and bear in mind that 
H(k) is easily convertible to Rij(k) using (24 1. 

A more detailed analysis of all these relations can be found in Moffatt ( 1978] >. 

For the magnetic energy density in ID Fourier space, a broken power-law is assumed in the following in our examples by 
adopting 



(25) 



usually with /3 = 2 and Icq = 1 if not stated otherwise, but with different spectral indices a. The low-A; asymptotic eg ~ k 2 
corresponds to a white noise spectrum without correlations on scales larger than 1/fcn. For large k, we find e B oc k~ a , eventually 
becoming a Kolmogorov-spectrum for a — 5/3. 

/v 2 2 

When ever necessary, we can always model the helicity power spectrum as H(k) = = |j-/z(£)e B (fc), where h(k) is a 

function between - 1 and 1 . This can be seen from 



My(k) 



e B (k) 
k 2 



k(k j 



^ 2 ) - in e ijm h(k) — 



(26) 



and the fact that the matrix Ay must be positive definite. We adopt k = ke x without loss of generality and find the characteristic 
polynomial of Ay to be 



(1 -Af-h 2 = 
— > 1 + h = A > 
— » \h\ < 1 



(27) 



This yields that h e [-1, 1]. 



2.3. General framework for all calculations 

The correlation functions of our observables can be calculated in a general and consistent way which we want to present now. 

Before we start, some general remarks about the mathematics are in place. Throughout this study an expression such as j'B 
relates to a multidimensional scalar product: 



>b = y f 



dx" J*(x)Bi(x). (28) 

i 

This definition includes a discrete summation over indices as well as a continuous integral over position space. The symmetric 
properties of matrix objects defined over a space with a scalar product ( |28] l reflect the appearance of discrete summation and 
continous integration. Therefore, a matrix element My(x, y) is called symmetric (or hermitian for complex quantities), if it is 
symmetric under a transposition of its indices and under an interchange of its vectors r: 

Mj j (x,y) = M ji (y,x). (29) 

Thus, a symmetrised element is expressed as 

My, sym (r) = i(My(r) + M^-r)), (30) 

where r = y - x. In the case where a matrix element is only symmetrised for index transposition, we call it index-symmetric: 

M, 7 , isym (x, y) = ^(My(r) + M^r)). (31) 

This distinction between symmetric and index-symmetric is important, because it takes care of subtleties that could easily generate 
confusion. We like to emphasize the difference between both symmetry operations, when applied to the magnetic correlation tensor 



fiSh. The tensor contains an intrinsic symmetric and an intrinsic antisymmetric element. Regarding ( 30 1, the intrinsic antisymmetric 



part is preserved, whereas regarding ( 3 1 1 it is not. This is of paramount relevance, since information on the helical power spectrum 



is only preserved, if the intrinsic antisymmetric parts do not cancel out during calculations. 
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Furthermore, we like to introduce the functional derivative, which is the natural generalisation of a derivative to function vector 
spaces. Its precise definition is (see Peskin & Schroeder|1995 i: 

3j(y) = S\x - y)6 ij 
<5 



£I/(x) 

6 1 
dJi(x) 2 

= J> 

1 



jtMJ= sTix) f dy3 f ^^^Maiy^W) 



i i 

-M, 7 (x,y')J/(y ) + ^Jk(y )M ki (x,y') 
1 



2 -(J r M) ! -(x)+-(MJ),-(x) 

For convenience and to avoid confusion with the delta function, we sometimes adopt easier notations: 

° = 3 y ,(x) = 3,-(x) 



c5J,(x) 



(32) 



Now we proceed, presenting the framework of the calculations. The general evaluation of the expectation value of a function X 
of observables for Gaussian magnetic field statistics with covariance matrix M and its determinant \M\ is conducted as follows^] 

<X(k ± ,kl,...)> B = -=L= f 2)BX(k ± ,kl,...)exp[-iB t M- I B] 
y/\2nM\ J 2 

exp[--B t M" 1 B] exp[/(kj_x ± + k' ± x' ± + ...)] 



),-.) 



- vpbi I DB fa - fa f dz ■ f dz "' F ^ (x) ^ (x '- ) ) 

exp[--B f M _1 B + JB] exp[/(k ± x ± + k^ + ...)] 
= J dx± ... J dx'" J dz ... J dz'" exp[/(kj_Xj_ + k^x^ + ...)] 

F (<9 y ,.(x), ^.(x'...)) | J=0 exp[i J'MJ] 
= J" dx ■■■ j " dx'" exp[/(k ± x ± + k^x^ + ...)] 

F (dj,(x), <9y,(x'...)) | J=0 exp[i J'MJ] 



J=o 



(33) 



(34) 



At this point, it is necessary to consider the actual form of F. The calculations vary enormously as complexity of the correlation 
functions differs. From now on we identify the covariance matrix of the Gaussian distribution with the magnetic correlation tensor. 
This identity is proven in Appendix |C| In ( |33] l we introduced a generating functional J and completed the square of the exponential. 
By integrating out the shifted Gaussian part, which depends on B, we are left with the part solely dependent on J: 

- 1 f £)Bexp[--B 1 "]Vr 1 B + JB] 
Vj2iM| J 2 

i r l .,. , l + l + 

= £>Bexp[--B , ]Vr 1 B+ -JB + -B'j] 

y/\2^M\J 2 2 2 

1 

~ Vi2^M| 



exp^j'MJ] 



f £)Bexp[--B t M" 1 B + -j'MM^'B + -B'M^'MJ] 

\nM\J 2 2 2 

= f ©Bexpt-^B-Mjy'M-^B-MJ)] exp^j'MJ] 

'.7iM\ J 2 2 



(35) 



In ( 34 1, we finally made the assumption that the observed space is sufficiently extended, so that we can neglect the finiteness of the 
integrals over the lines of sight J dz and treat them as if they were infinite. 



1 We denote with X either /, P or (f> or combinations thereof up to fourth order in the magnetic field in Fourier space. Therefore X contains up 
to 4 Fourier vectors within the observed plane labeled with primes. We denote with F the real space source function of X which depends directly 
on the local components of the magnetic field S, so that X = Jdz... f dz"'F(B, B', B", B'"). 
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ionai 



Introducing the generating functional J and thereby changing the fields B(x) to the functional derivative (9j ; (x) = gj:(x) in i 
provides a powerful method to calculate the integral DB over all possible magnetic field configurations as an infinite-dimensional 
path integral. With the definition of the functional derivative at hand < [32] >, we can discuss its actual evaluation. Since all uneven 
products of functional derivatives in p4| give zero because they also leave expressions with J which have been "taken down" from 
the exponential function during the differentiation, we are left with just two general types of possible combinations: Two or four 
derivatives. 

In addition to the generating functional technique familiar from quantum field theory, we also apply the renowned Wick theorem 
(see for example [Peskin & Schroeder|1 995 ) to evaluate the remaining derivatives in an elegant, quick and safe manner, rather than 
calculating them by brute force. The Wick theorem can be used under conditions which will become clear if one looks at how the 
underlying differentiation works. Firstly, we need the covariance matrix to be symmetric or hermitian. This, as already mentioned, 
means 



M !7 (x,y) = Mj/x,y) = M^(y,x) 



(36) 



which is fulfilled by (|7| and (Rjl. However, since there will be one case where the covariance matrix is not symmetric or hermitian 
(see Chapter[3]) and since we have to take thorough care of the exact order of the vectors, discussed in more detail below, the Wick 
theorem is expressed in a form that takes care of these subtleties: 



<9,(a)<9j(b) exp 
1 



- 2 rw 



.1=0 



2 (J , M)/b)+-(MJ) > (b) 









exp 




J=o 



(37) 



and 



3,(a)3,(b)^(c)3 / (d) exp 
1 



rJMJ 



.1=0 



2 (J , M) / (d)+-(MJ) / (d) 



<9;(a)d 7 (b)d*(c) 
3,(a)5 / (b)[(^M )t/ (c,d) + ^M tt (d,c)) 







exp 





.1=0 



(i(MJ),(d) + i(J 1 M) / (d))(^(J t M),(c) + \{MJ) k {cj) 



"2 
= 3/(a) 







exp 





.)=() 



(^M M (c,d) + iM tt (d,c))(i(J t M) ; (b) + ^(MJ)/b)) 



+ (^M y (d,b) + ^M^b,d))(VM),(c) + ^(MJ)/c)) 



^2 

+ (htfMMd) + ^(MJ) / (d))(^M y (c,b) + ^M^(b.c) + 



2 

II 
2 







exp 





= (^M H (c,d) + jM tt (d,c))(^(a,b) + ^M j7 (b,a)) 
+ (^M y (d,b) + ^Mj,(b,d))(^M ki (c,a) + l -M ik (*,cj) 
+ (^M, 7 (a,d) + ^M u (d,si))(^M k j(c,b) + ^M Jk (b,c)). 



J=o 



(38) 



Now all derivatives up to fourth order can be calculated just by simply inserting the current case. For example 
<9j(x)3j(x')|j=o exp ^j'MjJ can be read off from (38 i by inserting i = j = k = I = l,a = b = x and c = d = x' : 



di(x)dl(x') exp 



-JMJ 



= M 2 u (0) + 2M 2 u (r) 



(39) 



.1=0 



As we can see from ( |38| ), all antisymmetric parts cancel out during the differentiation. Thus although M was not explicitly restricted 
to be symmetric, only the symmetric elements of the magnetic correlation tensor 



M Usym (r)= Um u (t) + Mjt-T)) 



(40) 
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remain in the end. However, it is important to understand that these symmetric elements actually preserve the intrinsic antisymmetric 
parts that constitute the magnetic correlation tensor. As mentioned above, this is because we take into account the inversion of the 
vector r when transposing the tensor elements. A look at (|7j) reveals, that the minus sign of the Levi-Civita-tensor e ijm we encounter 
under interchanged indices is exactly cancelled by the minus sign occuring due to inversion of the vector r: 



M y (r) = M N (r)6 u + (M L (r) - M w (r)W + M H {r)e ijm r, 



(41) 



This means that although the tensor (41 



is symmetric in the general way defined in (30 1, it is not index-symmetric due to its 



individual antisymmetric constituents. The consequence is that if we carry out the derivatives using the Wick theorem we not 
only have to take care of the right combination of indices but also of the corresponding vectors r or -r and, in the end only the 
symmetric parts of M, as defined in ( [30] >, appear. And further, this does mean that if we encounter index-symmetric expressions 
such as Mjj(r) = (Mij(r) + M ; ,(r))/2, the intrinsic antisymmetric part related to the helical power spectrum in (41 1 is lost during the 
differentiation. A careful look at p8j ) reveals that with the right combination for i, j, k, I and a, b, c, d and a sum of terms as in (38 1, 
it is possible to get such combinations. For example 



(«9?(x)^(x') + ^(x)d?(x'))exp 



1 



rJMJ 



| J= o = 2M„(0)M 22 (0) 



+ 2(M|[(x-x') + M 2 2 (x-x')) = 2M„(0)M 22 (0) + 4M 2 \ lsym (x-x') 



(42) 



3. The polarisation 2-point function (P(k ± ) P*(k' ± )) B 

In this section the polarisation 2-point function (P(k ± ) ■ P*(k' A ))B is calculated to serve us as an example for the general calculation 
to obtain the other correlation functions of our observables. Since the steps are similar for all correlation functions and differ only 
in complexity, we intend to present them in detail only for a single case here and just list the other calculations in the Appendix [A] 
As (P{k x ) ■ P*(k' 1 ))B is of fourth order in the magnetic field and, in addition, P has a rather complex dependence on B, it is 

convenient to introduce a compact notation for P(x ± ) = L dz(B\(x) + /B 2 (x)) 2 in order to clarify the calculation as far as possible. 

Defining B ± = ^(B\ + iBq) allows the expression P(x x ) = dz 2B\(x) and P*(x x ) = j£ dz 2B 2 (x). Thus, a change of basis of 

B is introduced, mapping B = (Bi, B 2 , B3) — > B = {B+, B , B3). We can then work effectively with 2Z? 2 (x) • 2B 2 (x') instead of 
(Z?i(x) + iB 2 (x)) 2 ■ (Bj(x') - iB 2 (x')) 2 . With regard to the correlation function, this results in: 



<P(k_0 • P*(k'J) B = 4 j dx x j dx' ± Jdzj dz'exp[i(k x x x -k' x x' x 
d 2 J+ (x)d 2 j (x') exp 



)] 



rJ MJ 



(43) 



Regarding the differentiation with respect to J ± we need to establish a relation between J1/2 and J ± . The basis transformation should 
preserve all scalar products, therefore, we have j{B + + jlB- + J^B^ — J\B\ + /jfi 2 + jtBj,. Using this we establish the required 



relations: 



J^B+ + J B — J^B\ + J1B1 
1 



V2 

from which we obtain 
1 



(/| + JL) and /' 



V2 



(4 - /I) 



and also 7. = 



Thus, the transformation matrices J = OJ and J 



O 



J± = 
(J- + /+) 



1 



V2 
O 1 J ar e 



V2 
as well as 



h = 



V2 



(J- ~ J + ). 



(44) 

(45) 
(46) 



/ 1 


1 


0^ 




f 1 


i 


(T 


V2 






VI 


VI 


/ 







, and O 1 = 


1 


I 





V2 


VI 


VI 


VI 


I 





1, 




lo 





1, 



We now need to express the argument J' MJ of the exponential in (43 1 in terms of the transformed quantities: 

JMJ = JOMOJ = JMJ with M = OMO 
Some elements of M that will soon become important are: 

M ++ (r) = l - (M„(r) + M 22 (r) - iM 12 (r) + iM 21 (r)) 



(47) 
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M__(r) = l - (M„(r) + M 22 (r) + iM n (r) - iM 21 (r)) 
M + .(r) = i (M n (r) - M 22 (r) + zM 12 (r) + iM 21 (r)) 
M_ + (r) = I (M n (r) - M 22 (r) - /M I2 (r) - /M 21 (r)) 



The entire matrix then reads: 
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(48) 
(49) 
(50) 



M = 



M ++ (r) 
M_ + (r) 



M + _(r) 



■4=(Mi3 + M 23 ) 



M__(r) ^(M 13 -M 23 ) 



V V2 



^(Af 3 i - ;'M 32 ) ^(M 3 i + z'M 32 ) M 33 



Returning to the correlation function (43 1, we can now carry out the functional derivatives. Since J is a complex quantity, we now 
have to concern ourselves with the complex conjugation implied in the t operation which affects J = (/* , J*_, JV) = (J-, J+, 7 3 ). 
Using (38 i we find : 



d J+ (x) d J+ (x) djjx') djjx') expt^j'MJ] | J=0 = 2(M„(r) + M 22 (r)) 2 



(51) 



Inserting this into the overall equation for the correlation function (43 i results in: 



(P(k ± ) ■ P*(k'J) = 8 Jdx 3 J dx 3 '(M u (r) + M 22 (r)) 2 exp[/(k ± x ± - k' ± x'J] 

= 8 j dx 3 j dr 3 ■ (M 2 ,(r) + M 2 2 (r) + 2M 22 (r)M n (r)) exp[/x ± (k ± - k'J] 
exp[-ir ± k'J 

= 8 j dx 3 j dr 3 J J" -^.(M 11 (q)M 11 (q') + M 22 (q)M 22 (q') 

+ M 22 (q)Mn(q') + M 22 (q')$n(q)) exp[;x ± (k ± - k^)] exp[-ir(q + q')] 
exp[-;r ± k^J 

= 8(2;r) 2 5 2 (k ± - kl) j dz J dr 3 J ^ j ^ ■ (Mn(q)Mn(q') 

+ M 22 (q)M 22 (q') + M 22 (q)M n (q') + M 22 (q')M u (q)) exp[-/r ± (q ± +q' ± + k'J] 
exp[-ir z (q z + q'J\ 

= ^ 2 ( k ± ~ K)Lz J dq 3 J dq' 3 ■ 5 2 (q ± + ql + k'J6(q z + <?;)(M„(q)M 11 (q') 

+ M 22 (q)M 22 (q') + M 22 (q)M n (q') + M 22 (q')Mn(q)) 
a=(_<u = kl '" ?z) A<5 2 (k ± - k'JL z j dq 3 ■ (M„(q)M 1I (a) + M 22 (q)M 22 (a) 

+ M 22 (q)M n (a) + M 22 (a)Mn(q)) 

'^JV ^ [('-!)(' -I) 



In 



(2nfd 2 (k ± - k' 



S(2x) 5 d 2 (k x - k'JL z j 



dq 



e B (q)e B (a) 

q 2 a 2 



a- 

2 
a l 



(52) 



This integral can be further simplified. To do this, we transform it into spherical coordinates and perform the subintegral over ip 
analytically. The remaining 2-dimensional integral can be done numerically without problems. We choose the axes to be selected 
so that the angle 9 is between the x-axis and vector q ± while the angle tp rotates around the jc-axis. Without loss of generality, we 
choose k ± = k±e x to ensure that the angle between q ± and k ± coincides with 9. Thus the transformation is 



q x — q cos 9 

q y — q sin 9 sin <p 

q z - q sin 9 cos <p. 



(53) 
(54) 
(55) 
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Which implies for the required quantities: 



q\-q\ + q 2 - q 2 sin 2 6 sin 2 tp + cos 2 6 = q 2 (\ - sin 2 6 cos 2 tp) 



q 2 ± + k]_ + 2 q ± ■ - q 2 ± + kj_ + 2qkj_ cos 8 



a = (q ± + k ± ) + q\ = q\ + k]_ + q\ + 2 q ± ■ k ± 
= g 2 + fc 2 + 2gA: ± cos 0. 



(56) 
(57) 

(58) 



The integral (52 1 is then transformed as follows: 



K = ftf. ^^[(2-i)(2-4)] 



q 2 a 2 



dcos 6 



9 

e B (q)£B(a) 
q 2 a 2 



f 

Jo 



dtp [1 + sin 2 cos 2 y>] 



g 2 (l - sin 2 flcos 2 tp) + k]_ + 2qk ± cos 8 



I dq L' 



q 2 +k\+ 2qk ± cos 6 



dcos 9 



e B {q)e B {a) 



q 2 a 2 



2 + sin 2 6 



3 , „\/ ^sin 2 ^ 



+ (l + -sin 2 0)( 



<7 2 + fc 2 + 2gfc ± cos 0' 



(59) 



We integrate this numerically for values of k ± between k = 1CT 2 and 10 3 . We vary the spectral index a between a = | and the 
Kolmogorov-type spectrum a = 5/3. The results can be seen in Fig. |2](left). The slope of the declining section is not equal to a but 
depends on it. It is referred to as the polarisation spectrum slope a*. If one plots the energy spectrum slope a against the polarisation 
spectrum slope a*, one can see that there are two different regimes, for roughly a > 1 and a < 1 (see Fig. |2](right) ). 

It is now necessary to understand the approximate behavior of our findings in Fig. [2] (left) as well as the significance of the 
different regimes seen in Fig. [2] (right). As long as a is large enough the term e(g) only contributes to the integral within a sphere 
with radius qo in g-space, because there it is mainly constant: ^ w const. Beyond this sphere it is strongly suppressed by the 
(/""-dependence. In this case the slope is determined by the second term e(|q + k^J) in (59 1 which also contributes only within a 
sphere determined by q around the point k' ± . Inside this sphere, e is roughly eik'^) because q is small in comparison to k' ± since 
we are looking at the case k' ± > qo (otherwise we get a constant behaviour of the total integral as can be seen in the plots). So 
we get approximately / oc e(Q)e(k' ± ) which leads to log[7] * -(a + 2)log[A^] + const. This is confirmed by the approximations 
plotted in Fig. |2](left) and Fig. |2](right). In Fig. |2](left) we plotted a rough estimate for the integral, where we have just integrated 
I oc €(Q)e(k' A ) inside the go-sphere. In order to match the original integral better, it had to be shifted by a factor of 1 .6, which is 
perfectly reasonable considering that the simple sphere is just an approximation for a more complex structure. In Fig. [2] (right) we 
see that -a - 2 is indeed a good approximation for a* in the high-cr regime. The regime where a < 1 is not really of physical 
interest because all energy spectra with a < 1 would lead to the unphysical situation of infinite energies on the smallest scales as 
the spectrum complies with k~ a . Therefore we are not interested in the exact behaviour of the integral values below a - I. 




-2 -1 1 2 3 ' 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



Figure 2. Left: Numerically evaluated integral values of {P(k ± )-P(k' ± )) B in a log-log diagram with an approximation for the a = 5/3 
case. Right: Plot of the energy spectrum slopes a against the polarisation spectrum slopes a* in (f(k ± ) ■ Pik'^g- 
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4. Other correlation functions 

We now provide the results for the other correlation functions of our observables without repeating the details of the calculations. 
All functions of first or third order in B are omitted because they are obviously zero due to the uneven number of fields appearing in 
them. These are (<p(k ± ))ft, (<f>(k±)(f)(k' ± )(p(k"))B, Wkj_)0(k^))B an d (P(^-±)<P(^'j_))b- As stated before, the calculational steps to gain 
these expressions are fairly similar to the case of (P(k ± ) • P*(k' 1 ))B- As a matter of fact, most of them are even easier to obtain. More 
complex numerical integrations are only needed in the cases of (7(k ± )/(k^))B and (7(k ± )P(k^))B. In the following, the vectors r or 
r' shall always denote a combination such as x' - x, to be defined for each correlation function in Appendix[A] By T>{a, b) we denote 
the Beta-function. Furthermore u, v and w are defined as w = (k"', 0), u = (k", 0), v = (k' x , 0) and a = (— q x - k' ± , -q z ). 
Here are the results: 



</(kJ> B 

</ J (kJ) B 
<^(k ± )^(kl)> B 

</(kJ/(k'J> B 



2(2n) 2 5 2 (k ± )L z M N (0) 



<P(kJP*(kl)> B 



<0(k ± )^(kl)0(k^)0(k^')>B 



</(kJP(kl)) B 



</(kJ0(kl)0(k^)> B 



<P(k ± )0(kl)^(k^)> B 



= 128^ 2 (k^(| + i 
= 

= (2n) 2 6 2 (k ± +k'JL z M 33 (k' ± ,0) 
= 32n 5 6 2 (k 1 _ +k'JL z £(v)/v 2 



(2n) 2 6 2 (k 1 _ + k'JL z J ^(2(M ll (q)M„(a) 
+ M 22 (q)M 22 (a)) + 4M 2Usym (q)M 21 

,isym (a)) 



8n 2 6 2 (k 



x + k'JL z J 



dq 



e(q)e(a) 



q 2 a 2 



('-!)(' -5) 



♦(•-PX'-SMW) 

= |^ 2 ( k ^ " K)h J" dq 3 

(M n (q)M„(a) + M 22 (q)M 22 (a) 
+ M 22 (q)Mn(a) + M 22 (a)M n (q)) 

£B(q)£B(a) 



= 8(2;r)V(k 



i± -k'JL z j 



dq 



q 2 a 2 



= (27r) 4 [c5 2 (k" + k'")6 2 (k ± + k'JM N (w)M N (y) 
+ 5 2 (k'" + k'J6 2 (k" + k ± )M N (w)M N (u) 
+ 6 2 (k'[ + k'J6 2 (k'" + kJM N (w)M N (u)] 

M N (q)M N (a) 



= (27T) 2 6 2 (k ± +k'JL z f -p- 
J (2ny 



q 2 a 2 



[2{q\ + ql)(a 2 2 + a\) + 2(q 2 + q\)(a\ + a 2 )] 

= L 2 (27r) 4 6 2 (k' ± + k")5 2 (k ± )M iv ( M )2M A ,(0) 
- 2L,(27r) 2 5 2 (k ± + k' ± + k'l)H(u)H{v)luv 
(miVi + U 2 V 2 ) 

= 2L z (2n) 2 S 2 (k ± +k' ± + k") 
H(u)H{v) I uv ((uiVi - m 2 v 2 ) + i(u\\'2 + «2 v i)) 



(60) 
(61) 

(62) 



(63) 



(64) 



(65) 



(66) 



(67) 



(68) 



Many of the results are, as expected, providing no surprises. The mean total intensity ( |60j > is given, in principle, by the energy density 
of the magnetic field whereas the mea n polarized intensitv $6l\ should b e zero due to the is otropy of the problem. The correlation 
function (62 1 was already evaluated by EnBlin & Vogt| ( [2003| ~and also by Cho & Ryu] ( |2009| , it depends on the k z — plane of M zz . 
The fourth order quantities ( 63 i, ( 64 1 and ( 65 i are more complex, but nevertheless, only correlated combinations of the second order 
quantities. They can be used to monitor the validity of the assumption of Gaussianity and isotropy. Non-Gaussianity or anisotropy 
in magnetic field statistics would lead to a deviation from this form, which by comparison to ((p(kj_)(p(k' ± ))s can be detected. 

By far the most interesting results are of course ( 67 1 and ( 68 1 as they contain a direct dependence on the helical power spectrum 
H(k). Both are plotted in Fig. ^ at the same time, but only for the case where k ± = 0. We also restrict us to k = k e v for (68 1 
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Figure 3. Helical spectra of (/(kj_)0(k^)0(k"))B and (P(k ± )0(k^)0(k"))B for a = 1 and a = 5/3, k x = and different values of 
h(k) and under further assumptions for which both functions take on the same analytical form. 



without loss of generality. This reduces both correlation functions to 2(2n) ft (k). Actually plotted is only H (k). To represent the 
helical spectrum ft 2 (k) graphically, we factor out the energy spectrum esik) thus leaving us a function h(k), that parameterises the 
plots (see Sec. (2.2 1). 

The exact results (67 1 and (68 i can be understood from the physical point of view. In both correlation functions we have a part 
connected to the magnetic field component that lies in the surface perpendicular to the line of sight /(k ± ) and P(kjJ due to the 
polarisation properties of synchrotron emission. There is also a part which depends on the line-of-sight component <p{k ± )<p(k' ± ) due 
to Faraday rotation. If we consider the corresponding delta functions, we see that in cases, where we have 6 2 {k' ± +k")(5 2 (k ± ) the two 
parts are somewhat uncorrelated as the related term is more or less the product of (7(k ± ))B and (</>(k ± )0(k^))B. This is reflected in 
the absence of any helicity dependence. In contrast, if we examine the parts with 6 2 (k ± + k' ± + k") there is a type of mixing between 
the different observables due to the mutual dependence of the three vectors k ± , k^_, k" through the delta function. Accordingly, there 
is a dependence on helicity. More about the physical interpretation of these results can be found in the following section. 

Please note that information on the overall sign of the helicity cannot be obtained using this method becaus H{k) only appears 
quadratically. 



5. The LITMUS test 



For the correlation function (68 I, a strikingly intuitive picture can be found to explain the result. Let us take a look at Fig. Q where 
we imagine the line of sight to be directly aligned with the axis of a magnetic helix. In a combined polarisation and Faraday depth 
map, we should see a central region with nonzero faraday depth (p and around it a radial polarisation pattern. These are correlated 



structures, that make (68 i nonzero for helical magnetic fields. They would vanish with the field becoming nonhelical. However, we 
clearly have two possibilities for the direction of the magnetic field going around and therefore could get positive or negative <p 
respectively. Thus, these correlated structures can only be seen in (P(k ± )0(k^)0(k"))B and not in (P(k ± )0(k^))B, where the single 
dependence on (f> would induce the positive and negative parts to cancel out over averaging. This is confirmed, as (P(k ± )0(k^))B 
becomes zero due to the odd number of functional derivatives. 

Guided by this picture, the LITMUS test (Local Inference Test for Magnetic fields which Uncovers heliceS) was developed, a 
small and simple test that could be easily used to probe real data for helicity. 

The LITMUS test was not constructed to produce quantitative measurements of the helicity spectra, but to provide a fast and 
qualitative test for the presence of helicity. 



Henrik Junklewitz & Torsten EnBlin: Imprints of magnetic and helicity spectra 



13 




Helical magnetic fieldlines 
in two different direction 



Figure 4. Schematic picture of correlated structures in combined polarisation and RM maps that can give rise to a non-zero corre- 
lation function <P(kj_)0(k^)0(k")) B . 



5.1. The basic idea 

We take a closer look at one of the aforementioned patterns of polarisation and Faraday rotation (see Fig. HJ. In the ideal, helical 
case, the gradient G = of (f> points either to the center or out of the center of the region with Faraday rotation and should therefore 
be perfectly aligned with the polarisation. The polarised intensity P is a complex number representing a spin 2 field. To compare P 
with G, we just transform G from a two-dimensional vector into a complex number Q in the same representation of "directionless 
vectors": 

Q = |G| 2 exp[2i'a] with a = arctan -p. (69) 

G x 

By doing so, we loose the information on the direction in which the gradient G is pointing. The quadratic dependence of @ on |G| 
accounts for a normalisation that will become clear in the following. 

We now want to construct a test using Q and P that is sensitive to the presence of magnetic helicity but can be performed easily 
on a real dataset. Our previous considerations suggest to simply use the fact that Q and P should be parallel for helical fields. The 
test then just consists of multiplying Q with P* for every pixel of a given map of (p 2 and P. This complex scalar product produces 
different results for different orientations of Q and P in the complex plane. If the gradient and the polarisation are parallel (% and a 
differ by a multiple of n) the result is real and positive. If they are perpendicular (x and a differ by an odd multiple of n/2) the result 
is real and negative. For any orientations in between, the result will take on complex values. 

We now can state the LITMUS test. In the presence of helical fields, the average of the scalar product Q*P over all pixels of a 
f^T-map should have a real value significantly larger than 0. Whereas in the case of non-helical fields, the alignment of Q and P 
should be changing randomly from pixel to pixel so that we would expect the average over Q*P to be zero. In short mathematical 
notation this is stated as 



(e*p) 

(OTP) 



> and real, (70) 

helicity 

= 0, (71) 

no helicity 



where the ensemble average over helical or nonhelical fields is in practice replaced by an average over all pixels of a £?*P-map. 

It can be shown that our intuition was right and that this test actually depends directly on (P(k)0(k')0(k")) which was our intial 
starting point. We begin by writing out the full condition for the LITMUS test: 

Iff (x)P(x))) = | [(^) 2 + (^) 2 j exp [ - 2i arctan G y /G x ]P(x^. (72) 
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Using trigonometrical theorems it is easy to show that 

(p 2 K + G 2 ^ exp [ - 2/arctanG v /G A ] = (p x - iG^j . 
We now apply a Fourier transformation to our observables P(x) and <p(x) and rewrite the whole expression: 



(73) 



{g*(x)P(x)} = 



dd>(x) 84>(x) 
— i 



dx\ 
d d 

dx-r l dx~2 



8x2 



P(x) 



<p(x) 

.2 



d d 

dx-r 1 ^ 

a" 



x)\P(x) 



P(k)0(k')0(k")) 

= / ^l^-^Pi-^] exp[/x(k + k' + k")] 

/p(k)0(k')0(k") 



(74) 



We assumed again that the average over all pixels of a §*P-map is equivalent to an ensemble average over the magnetic field 
statistics. 

Now we insert our result (68 1 for (P(k)0(k')0(k")) and see the dependence of ^(x)P*(x)^ on the helical spectra H(k')H(k"): 

2" 



{g(x)nx)) = 2L z (2n) 2 J ^^(k + k'+k")^-^]^'-^'] 



exp [z'x(k + k' + k")]^'^" - k^k^) + i(k[k' 2 ' + fc^fc")] 



H(k')H(k") 
tk" 



= 2L z (2n) 2 f ^ J ^[^-^][^-^'][(r i ^-^ 2 ') + /(fe' 1 ^'+^' 1 ')] 



H(k')H(k") 
tk" 



(75) 



, r dk 2 ' r dk 2 " , ,, - , . „ 



H(k')H(k") 
k'k" 



= 2L 7 



f 

Jo 



dk k 2 H(k) 



= 2L,7T 4 



f 

Jo 



dk 



e H (k) 



(76) 



In (75 1 we assumed w. 1. of g. that x = and in the last line we used ( 17 1 to substitute e#. We see that ^(x)P*(x)^ is a direct 
and clear estimator which measures the square of the k-space integrated helicity spectrum, giving more weight to the large scales. 
Therefore, we conclude that, except maybe from pathological or fine tuned situations in which the k-weighted helicities at different 
k-scales cancel each other (since H(k) might change sign), we can expect the LITMUS test to be able to reveal the presence of 
magnetic helicity. This is demonstrated in numerical tests on simulated data conducted by Oppermann et al. (2010). This work also 
contains an application on real data of our own galaxy and a thorough analysis thereof]^] 

Finally, we can state that, in principle, a given dataset can be tested for helicity. If P and <p are available, a test using the conditions 
(70 1 and (71 1 can easily be implemented due to the usage of local quantities only, avoiding any data transformation to Fourier space 
and the complication finite window functions thereby would introduce. 



6. Conclusions 

We have shown how statistical properties of turbulent cosmic magnetic fields can imprint on the statistics of certain radio observ- 
ables. Our analysis involved the total intensity 7(k ± ) and the polarised intensity P(k ± ) coming from radio synchrotron emission out 
of a volume as well as the Farady depth 0(k ± ) of background sources seen through the same volume. The first two depend on B ± , 
the magnetic field component lying in the surface perpendicular to the line of sight. In contrast, 0(k_i_) depends on By, the component 



2 Unfortunately, it seems that a straightforward application of this test can be hampered by a too large spatial variance in the elctron density, 
which we assumed to be constant in this work. 
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parallel to the line of sight. Whenever this set of observables is available, we can examine all three components of the magnetic 
field. 

With regard to these observables, we evaluated a complete set of cross-correlation functions up to fourth order in the magnetic 
field and presented simple analytical equations in Fourier space depending on the field's energy spectra. We demonstrated that two 
correlation functions of our set, namely (P(kj_)0(k^)0(k"))B and (7(kj_)0(k^)0(k"))B, explicitly depend on the helical spectra of 
the turbulent field. The first one depends solely on the helical parts and becomes zero for non-helical fields. 

This finding offers a new way for measuring the helicity of magnetic fields and thereby for testing existing mean field dynamo 
theories involving helicity. Measuring these correlation functions in real data will permit the study of helicity spectra up to their 



overall sign. If (67 1 or (68 1 provide a non-zero result with statistical significance for k ± + , this is direct evidence for helicity in 
the magnetic field. 

Furthermore, we presented the LITMUS test, a simple procedure to be applied to data which probes for helicity. The LITMUS 
test is easy to apply since it can be fully computed in real space. It provides the square of the k-weighted k-space integrated helicity. 
First results of an application of the test to real data and to numerical tests using simulated helical and non-helical fields can be 



found in Oppermann et al. (2010). 



Our general formalism permits the construction of further tests which can probe helicity on invidual k-scales. However, before 
such tests can be applied to real data, suitable observational configurations and further theoretical development in order to alleviate 
the simplifications and assumptions made beforehand are required. 

Subsequent work should therefore follow two directions: To find observations that match our assumptions best and to extend 
our calculations to be able to cope with more complex observational situations. Advancements should include 

- more realistic non-Gaussian components of the magnetic field statistics 

- the removal of statistical homogeneity as an overall simplification 

- spatially varying electron densities instead of assuming them as constant (see Appendix [B| 

- the possibility for preciser values for the spectral index of the cosmic ray electron density than the choice p — 3 (see Appendix 


- calculations without the restrictions imposed by observed space being large leading to the approximation J dz ~ J m dz 
(see Sec. |2.3) > 

- the introduction of a window function formalism 

- and finally developing an approach for dealing with intrinsic Faraday rotation which modifies the polarisation at long wave- 
lengths, a topic which has been neglected here (see Appendix IB), 

Before becoming more deeply involved in discussions on possible advances in the future, we should first think about which observa- 
tions could be applicable to our approach in the present form. Observations required by our analysis have to come from a polarized 
radio-synchrotron source with background Faraday-rotation. A suitable target could probably be found in the interstellar medium 
(ISM) within our own galaxy, of which we have some established knowledge in respect of large scale fields. With regard to the latter 
point, it would be advantageous to choose a region in which the magnetic helicity flows are expected to be found in accordance with 
mean field theory. Furthermore, the polarisation data has to be taken at high frequencies to be Faraday rotation free. For that, the 
upcoming Planck polarisation data of our Galaxy will be ideal, since it is at short wavelength, has high resolution and accuracy and 
is also full sky. Compilations of RM measurements of background sources seen through our galaxy already exist ( Haverkorn |2007 
Brown et al.|2007)|Taylor et al.|2009 1 



Attractive extragalactic objects to be investigated for magnetic field statistics are the lobes of radio galaxies, whose intensity and 
polarisation statistics can be constructed. However, no Faraday rotation could be detected yet through their lobes. For the radio jets, 
this is different and there, helicity can be probed and is yet actually expected to be present (EnBlin 2003; Gabuzda 2005; Mahmud 
|& Gabuzda|2008] >. 

Galaxy clusters are probably not well suited for our approach, although they host large scale magnetic fields. The high degree 
of intrinsic Faraday rotation usually found there erases polarisation of the cluster radio halo emission at the synchrotron frequencies 
we observe them. 

However, handling intrinsic Faraday-rotation analytically is a considerable challenge. We would need to include an extra expo- 
nential factor for the rotation in P(x ± ): 

P(x x ) = J*-[Bi(x) + /B 2 (x)] 2 -exp[2^(x ± )i 2 ]. (77) 

As 0(x ± ) itself contains an integration over dz, the extra exponential factor couples all positions along the line of sight. Thus, 
the exponential becomes so complicated that it has to be approximated in a suitable way. This seems to spoil a purely analytical 
approach. Nevertheless, the inclusion of intrinsic Faraday rotation is one of the next important challenges as it would enable our 
approach to be applied to many sources excluded until now, such as galaxy clusters. 

To resume further our basic discussion, there is more to consider. From a technical aspect, the first problem to be tackled for 
tests, which aim to measure detailed magnetic energy and helicity spectra, is the inclusion of a realistic window function into the 
formalism. A window function is set by the observations but also incorporates variation of the signal due to changing relativistic 
electron density and magnetic field strength. It scales with the electron density n e (x) for 0(k ± ), with the cosmic ray electron density 
n cre (x) for 7(k ± ) and P(k ± ) and with the average magnetic field profile. The two former densities have to be taken from independent 
observations (e.g. free-free emission), the latter has to be guessed depending on the source and relying on prior knowledge. It is 
clear that introducing a window function will make calculation and integration more complex. In principle, however, there is no 
basic restriction that could prevent its implementation. Previous attempts to measure magnetic power spectra from observations of 
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magnetic fields heavily affected from window functions have already proven successful, e.g. En Blin & Vogt| ( |2003] l; |Vogt & En Blin 
(2005|l; |Kuchar & EnBlin| ( |2009) . 

The next aspect to be considered is whether to use Gaussian statistics. Real fields are probably non-Gaussian. Any attempt to 
model exactly the real situation has to include at least non-Gaussian deviations. This brings additional complexity into our approach 



as they could not be handled analytically anymore. In ( 34 1, we would not get rid of the path integral and would need to rely on 
approximations or perturbative approaches, such as in field theory. 

At the beginning, our primary goal was, firstly, to find simple analytical relations between the statistics of radio observables 
and magnetic fields and, secondly, to prove conceptually that with such an approach, it is in principle possible to extract informa- 
tion about the helical part of the magnetic field. To show this, we therefore choose to start with Gaussian statistics, the simplest 
configuration possible, which also would allow us to keep our calculations analytical while higher order statistics always could 
be incorporated later in a perturbative expansion around the Gaussian case. We do not claim that this assumption is sufficient to 
reproduce exact results in accordance with the high complexity of real nature. However, we do believe that while the grade of 
Gaussianity in our statistics might determine the strength with which the statistics of observables depends on magnetic properties 
like helicity, the assumption of Gaussianity is not essential for the dependence itself to occur. Thus we believe we can, in principle, 
decide whether data contains signatures of helical magnetic fields or not and that the simplification of using Gaussian fields in our 
calculations is sufficient to achieve the goal of showing how helicity and other magnetic properties can be detected. This of course 
needs to be shown e.g. using mock data from numerical simulations, but this is left for further studies. 

The reminder of our assumptions and simplifications only represent minor problems. Deviations from the cosmic ray electrons 
spectral index p = 3 could be included in the form of correction terms. To assume that the distance between observer and source 
is very large (and therefore the line of sight projection is parallel) is most of the time a fairly well approximation given the vast 
distances we encounter on cosmic scales, but it might break down if we analyse the large scale magnetic field directly in front of 
us. For the LITMUS test, the work of Opperm ann et al.| ( |2010| l shows that it is still applicable in such a situation. However the 
varying n e seems to be a more severe problem. The assumption of statistical homogeneity is widely used in the literature and proved 
appropriate for similar problems in the past. 

The outcome of our study will hopefully contribute to new findings on cosmic magnetic fields. It establishes a new, structured 
and definite way to measure magnetic helicity and it offers a new option to test cosmic dynamo theories. The presented LITMUS test 



might be able to easily probe data for helicity. However, first results by|Oppermann et aL (2010 1 show us that further development 



is necessary for an actual helicity-sensitive implementation given the complication in realistic observational situations, especially 
the varying electron density. Last but not least, we present a complete range of correlation functions of radio observables in elegant, 
simple forms, which are easy to evaluate and, in principle, can be compared to real data. In fact, some of the correlation functions 
may be of interest in themselves, disregarding the topic of helicity. 
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Appendix A: All correlation functions 

This Appendix contains a full list of all evaluated correlation functions and their derivations. This constitutes the core of our work 
but provides no further inside information on the concepts or helps understanding the findings. Therefore, we have collected the 
calculations in this Appendix. The details of the calculations are similar to the example of (P(k ± ) ■ P*(k' ± ))B, dealt with some depth 
in Sec. [3] and are only commented on if necessary. All correlation functions are given in Fourier space. The functions with an odd 
number of fields are equal to zero due to the analogous odd number of derivatives with respect to the generating functional J. The 
calculations involving P(k ± ) are given in a Faraday-free case, but set up in a way so that we could include Faraday rotational effects 
for further analysis later on. Refer to chapter[6]for more information on future perspectives. In the following, the vectors r or r' shall 
always denote a combination such as x' - x, to be defined for each correlation function in Appendix [A] In the following u, v and w 
are defined as w = (k"',0),u = (k",0), v = (k^,0) and a = (-q ± - k^_, -q z ). 

A.1. Calculation of {<f>(k ± )) B 




(A.l) 



A.2. Calculation of {<f>{k ± )(f>(k ± ))B 



<</>(k±)0(kl)>B = I dx 3 I dx' 3 exp[ik ± x ± + ik^x'J 3 3 (x) 3 3 (x') 




exp -J'MJ | J=0 



= (2n) 2 S 2 (k ± + k'JL : J dr 3 exp[ir JL k'J M 33 (r) 
= (2;r) 2 <5 2 (k ± + k'JL z J dr 3 J CK V [ir ± k'J exp[-/rq] M 33 (q) 

= (2n) 2 S 2 (k ± + k' ± )L z j dq 3 6 2 (k' ± - q ± ) 6{-q z ) M 33 (q) 

= (2n) 2 8 2 (k^ + k'JL z M 33 (k' ± ,0) 
= (27r) 2 5 2 (k ± + k'JL z M N {v) 

= 327!- 5 5 2 (k L +k' ± )L z ^ 




(A.2) 



v 



A.3. Calculation of {<f>(k ± )(f>(k' ± )(f>(k'l)) B 




3 3 (x)|j =0 3 3 (x')|j=o3 3 (x")exp -j'MJ | j=0 = 



(A.3) 



A.4. Calculation of {<f>(k ± ) t f>(k' ± )cf>(k'l)<p(k'")) B 



<0(k ± )0(kl)0(k^)0(k7)> B = n\ dx 3 dx' 3 dx' 
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exp[*k ± x ± + ik' ± x' ± + *X + iK'x^'][d 3 (x)d 3 (x')d 3 (x")d 3 (x'")] 
expjij'MjJ |j =0 

= n 4 j dx 3 j dx' 3 j dx" 3 j dx'" 3 exp[/k ± x ± + ik' ± x' ± + *X + ik'['x'['] 
[ M 3 3(x"'-x")M3 3 (x'-x) + M 33 (x"' - x')M 33 (x" - x) 

Parti Part2 
+ M33(X W -X)M33(X"-X') ] (A.4) 

Part3 

Part 1 gives with x'" - x" = r' and x' - x = r: 

Part 1 = n 4 J dx 3 J dx" 3 J ' dr 3 J dr' 3 exp[/x"(k" + k"')] exp[/x ± (k ± + k^)] 
exp[/r^k"']exp[/r ± k^]M33(r')M33(r) 

exp[ir ± (k' ± - qj] exp[-ir z q z ] exp[/r^(k"' - q^)] exp[-ir' z q' z ] 
M 33 (q')M 33 (q) 

= ^(27r) 4 5 2 (kl' + K')S 2 (k ± + kl)M 33 (k7, 0)M 33 (kl, 0). (A.5) 

Part 2 and Part 3 essentially provide the same if one just adopts the definition of r and r' as follows 
Part 2: x'" - x' = r' and x" - x = r, 
Part 3: x" - x = r' and x'" - x = r. 
The final result is then 

<<A(k ± )0(kl)0(k^)^(k7)) B = n 4 e (2n) 4 [d 2 (K + k"> 2 (kx + k^) 

M 3 3(k^',0)M 3 3(kl,0) + (5 2 (k7 + k'J6 2 (k'^ + k ± )M 33 (k'^' ,0)M 33 (k'^0) 
+ 6\k'l + k'J6 2 (k'l' + kJM 33 (k'[', 0)M 33 (k'[, 0)] 

= n 4 (2n) 4 [s 2 (k" + k'")6 2 (k ± + k' ± )M N (w)M N (v) + 5 2 (k"' + k' ± )6 2 (k" + kj 
M N (w)M N (u) + 5 2 (k'l + k'J6 2 (k"' + k ± )M N (w)M N (u)] (A.6) 

A.5. Calculation of (I(k ± )) B 

</(kJ> B = J" d /x 3 e X p[/k i x ± ](3 2 (x)+3 2 (x))exp ij f MJ | J=0 

= J dx 3 exp[i'k ± x ± ] (Mn(x,x) + M 22 (x,x)) 

= j dx 3 exp[/k L xJ(M 11 (0) + M 22 (0)) 

= 2M N (0) J"c/x 3 exp[/kj_Xj_] 
= 2(27r) 2 ( 5 2 (k ± )L,M A ,(0) 

= l2Sn 4 6\k ± )L z B^ + 1 -,^- 1 -) (A.7) 

In the final step, M N (0) was expressed in terms of a Beta-function S(a, b). This function is assigned a specific value for a given set 
of spectral indices a and p. 

A.6. Calculation of {I(k 1 _)I(k , J} B 

</(k J/(kl)> B = Jdx 3 J dx' 3 exp [i(k ± x ± + klxl)] (d 2 (x) + <9 2 (x)) 

(^(xO + ^Cx'^expfijtMjl |j=o 
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= jdx 3 j dx' 3 exp [i(k ± x ± + klxl)] (^(x)^(x') + <9?(x)<^(x') 



1 



rJMJ 



|J=0 



+ <^(x)^(x') + <9 2 (x)<5 2 (x'))exp 
J" dx 3 j dx' 3 exp [iOkjXj. + k' x x'J] ( (M„(0) + M 22 (0)) 2 

Parti 

+ 2(M 2 ,(r) + M 2 2 (r)) + 4M 2 ljSym (r) ) 



Part2 



Part l[^]is not dependent on r and can be applied directly without difficulties: 

j dx 3 j dx' 3 exp [i(k ± X ± + k x xl)] (Mn(0) + M 22 (0)) 2 

= (27r) 4 6 2 (k x )6 2 (k'JL z L' : (M u (0) + M 22 (0)) 2 
= (27r)V(kJ5 2 (k ± )L z L:4M 2 (0) 
In contrast, Part 2 needs a little more work using Fourier transformations: 

J" rfx 3 J dr 3 exp [/x ± (k ± + k'J] exp [ik> x ] (2(M 2 ,(r) + M 2 2 (r)) + 4M 2 2 l sym (r) 
= (27r) 2 5 2 (k ± + k'JL z J dr 3 exp [/k ± rj (2(M 2 ,(r) + M 2 2 (r)) + 4M 2 , sym (r)) 

= (27r) 2 5 2 (k ± + kl)L, fdr 3 J ^ J ^ exp [i\' ± r ± ] exp [fr(q + q')] 
(2(M„(q)M 11 (q') + M 22 (q)M 22 (q)') + 4M 21 , sym (q)M 21 , sym (q')) 

= (27r) 2 <5 2 (k ± + k'JL z Jdr 3 j^j^ exp \ir ± (q x + q' x + k'J] 
exp [zr z (& + ? ;)] (2(M 11 (q)M 11 (q') + M 22 (q)M 22 (q)') 
+ 4M 21iSym (q)M 21 , sym (q')) 

= (2*) 2 * 2 (k ± + k' x )L z f£f£ (2n) 3 6 2 (q ± + ql + k x ) 

<% z + ^)(2(M n (q)M n (q') + M 22 (q)M 22 (q)') + 4M 2I , sym (q)M 21 , sym (q')) 
= (27r) 2 5 2 (k ± + kl)L, J ^(2(M ll (q)M I1 (a) 
+ M 22 (q)M 22 (a)) + 4M 2 i >sym (q)M 2 i jSym (a)) 



= 8n 2 6 2 (k ± + k 



q 2 a 2 



jJ'MJ 



This integral has can be solved numerically for specific choices of the spectral indices a and yS. 

A. 7. Calculation of (P(k) J B 

(P(kJ) B = J dx 3 exp[/k ± xj (d\(x) - 8 2 (x) + 2/5i(x)3 2 (x)) |j exp 

= J dx 3 exp[/k ± xj Jm 2 ,(0) - M 2 2 (0) + 2/(^M 12 (0) + ^M 21 (0)) 
= 

The last step holds because of the relations M n (0) = M 22 (0) and M 2i (0) = M i2 (0) = 0. 
3 This is actually the same as the results obtained by computing </(k ± ))| 
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A.8. Calculation of {I(k ± )P(k ± )> B 

</(k ± )P(kl)) B = jdx 3 j dx' 3 exp [i(k ± x ± + k' ± x'J] (d 2 (x) + d 2 (x)) 

[d\(x') - d\{x' + 2/S 1 (x')d 2 (x')))exp ^MJ | J=0 

= j dx 3 j dx' 3 exp [i(k ± x ± + k^x'J^M^O) - M 2 22 (0) + 2M 2 n (r) + 2M 2 22 (r) 

+ 2M 2 l2 (r) - 2M\ x (f) + 4iM 2l (r)(M u (r) + M 22 (r))] 

= (2n) 2 6 2 (k 1 _ + k' ± )L z jdr 3 j ^ J ^ exp [ik' ± r ± ] exp [/r(q + q')] 

[ + 2M n (q)M n (q') + 2M 22 (q)M 22 (q') + 2M u (q)M u (q') - 2M 2l (q)M 2l (q') 
+ AiM 2l {q'){M n {q) + M 22 {q))\ 

= (2n) 2 6 2 (k ± + k'JL z f^f^ i^S 2 ( q± + ql + K) 

S(q z + q z )[2M u (q)M u (q') + 2M 22 (q)M 22 (q') + 2M l2 (q)M l2 {q') 
- 2M 21 (q)M 2l (q') + 4iM 2l (q')(M u (q) + M 22 (q))] 

= (2n) 2 6 2 (k 1 _ + k'JL z J ^-[2M n {q)M n {a) 

+ 2M 22 (q)M 22 (d) + 2M l2 (q)M 12 (o) - 2M 2l {q)M 2l {o) 
+ 4iM 2l (a)(M n (q) + M 22 (q))] 

= (2n) 2 S 2 (k ± + k'JL z J ^ f N( f 2 f (a) [2(q 2 2 + q\ ){ a 2 2 + a\) 



2(^i + q\)(a\ + a\) - 4i(q\ +q 2 2 + 2q 2 3 )a 2 a 3 ] + 2 



M N (q)H(a) 



q 2 a 



[iqiq 2 a 3 - 2{q\ + q\ + 2q\)a 3 } 



M N (a)H(q) 



q 2 a 



[icna 2 q 3 \ 



= (2n) 2 5 2 (k ±+ k'JL z f 
+ 2{q\ + q\){a\ + of)]] 



3 tM N {q)M N (a) 
) 3 ' 



q 2 a 2 



[2(q 2 + q 2 )(a 2 + a 2 ) 



(A. 12) 



The last step was possible using an asymmetric property of the integrals. All terms in which a 3 occurs become zero when integrated 
from -00 to 00. We can split up the integral into two parts covering the negative and positive regions and it can be seen that they will 
cancel each other. The remaining integral can be solved numerically for specific choices of the spectral indices a and /}. 



A.9. Calculation of (I(k ± )<f>(k' ± )) B 



(/(kj_)0(k^)) B = n e «o J dx 3 J dx' 3 exp[/kj_x ± + ik' ± x' ± ] 
(^(x)d 3 (x') + ^(x)d3(x'))|j =0 exp 



J MJ 



= 



A.10. Calculation of (P(k 1 _)4>(k'J) B 

(P(kj_)0(k^)) B = n e a Q J dx 3 J dx' 3 exp[;k ± x ± + ik' ± x' ± ] 

(d 2 (x)d } (x') - d 2 (x)d 3 (x') + 2/d 1 (x)d 2 (x)d 3 (x')) b=o 



exp -J'MJ 



(A. 13) 
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A.11. Calculation of (I(k A _)4>(k'J ( f>(k'l))B 

(/(kj_)0(k^)0(k")> B = J dx 3 J dx' 3 J dx" 3 exp[/kj_x ± + ik' ± x' ± + ?k"x"] 

(d 2 (x)d 3 (x')3 3 (x") + d 2 (x)d 3 (x')3 3 (x")) |j=o exp ^MJ 
= j dx 3 j dx' 3 j dx" 3 exp[/k ± x ± + ik' x x' x + ik"x"][ M u (Q)M 33 (x" - x') 



Parti 

+ 2 (M 31 (x" - x)M 31 (x' - x)) + M 22 (0)M 33 (x" - x') 



+ 2(M 32 (x"-x)M 32 (x'-x)) ] 

Pait4 

Part 1 and 3 can be calculated in the same way, as they only differ in the component of M chosen by the derivatives (i 

J dx 3 J dx' 3 J dx" 3 exp[/kj_x ± + ik' ± x' ± + /k"x"]M,/(0)M 33 (x" - x') 

x = x +r M,,(0) J dx 3 J dx' 3 J dr 3 M 33 (r) exp[/x^(k^ + k")] exp[/rj_k"] 
exp[/k ± x ± ] 

= M u (P)L 2 z (2n) 4 S 2 (k' x + K)6 2 (k ± ) J" d 3 rM 33 (r)exp[ir ± k' ± ] 

= M, 7 (0)L?(27r)V(kl+kl')<J 2 (kJ J d 3 r J J^LlCf 33 (q) 

exp[/r ± (k" - q J] sxp[ir z q z ] 
= M, 7 (0)L 2 (27r)V(kl + kl')<J 2 (k ± ) J dq 3 M 33 (q)6 2 (K ~ lJS(-q z ) 
= M u (0)L 2 (2n) 4 6 2 (k' ± + kl> 2 (k ± )M 33 (k^, 0) 
Equally, Part 2 and Part 4 can be solved on the same basis (i = 1,2): 

2 j dx 3 j dx' 3 j dx" 3 exp[/k ± x ± + ik' ± x' ± + ik'[x'[]M 3i (x" - x) 
M 3i (x' - x) 

= 2 J" dx 3 J" dr 3 J dr' 3 exp[/x ± (k ± + k^ + k")] exp[/r^k"] exp[/r ± k^] 
M 3i (r')M 3i (r) 

= 2U2n) 2 6 2 (k ± +k' ± + k'i) j dr 3 j dr' 3 j J exp^r^ - q^)] 

exp[-ir z q' z ] exp[;r ± (k^ - q ± )] exp[-ir z q z ]M 3i (q'))M 3i (q) 
= 2L z (2n) 2 6 2 (k ± +k' ± + k")M 3i (k", 0)M 3i (k' x , 0) 

This gives as final result: 

<7(k x )0(kl)0(k^)>B = L 2 (2;r)V(kl + k^)<5 2 (k JM 33 (k^, 0) 

[M„(0) + M 22 (0)] + 2n 2 L z (2n) 2 6 2 (k ± +k' ± + k^)[M 31 (k^,0)M 31 (kl,0) 
+ M 32 (k^,0)M 32 (kl,0)] 
= L 2 (2n) 4 6 2 (k' x + k")6 2 (k x )M N (u)2M N (0) 

- 2n 2 L,(27i) 2 5 2 (k 1 _ +k'+ k'') —^—^- (uiV } + u 2 v 2 ) 

uv v ' 

A discussion and a plot of this important result can be found in Sec. [4] 



22 Henrik Junklewitz & Torsten EnBlin: Imprints of magnetic and helicity spectra 

A.12. Calculation of (P{k x )^{k x )^{kl)) B 



<P(k ± )0(k^)0(k")> B = j dx 3 j dx' 3 j dx" 3 exp[/k ± x ± + ik' ± x' ± + zk^x"] 
(<9 2 (x)c) 3 (x')3 3 (x") - d 2 (x)d 3 (x')d 3 (x") + 2id l (x)d 2 (x)d 3 (x')d 3 (x")) | J=0 



exp 



= J dx 3 J dx' 3 J dx" 3 exp[/kj_x ± + ik' ± x' ± + ik"x"] 

[m 33 (x" - x')(M„(0) + M 22 (0)) + 2M 3I (x" - x)M 3I (x' - x) 
+ 2M 32 (x" - x)M 32 (x' - x) + 2iM 32 (x" - x)M 3 j(x' - x) 

+ 2iM 32 (x' - x)M 31 (x" - x)] (A. 19) 



This can be evaluated in exactly the same way as (A. 18 i because all the terms have the same basic structure. This provides: 

<P(k ± )0(kl)^)) B = L 2 (27r)V(kl +k^ 2 (k ± )M 33 (k^0)[M 11 (0) 
+ M 22 (0)] + 2n 2 e L z (2n) 2 S 2 (k ± + k' ± + k") 

[M 3 i(k^0)M 3I (kl,0) +M 3 2(k",0)M 32 (kl,0) + /M 32 (k^,0)M 31 (kl,0) 
+ /M 31 (k^0)M 32 (kl,0)] 

= +2n 2 „L z (2n) 2 6 2 (k ± + k' + k") — — (ttiVi - u 2 v 2 + i(u x v 2 + m 2 vi)) (A.20) 

uv v ' 

A discussion and a plot of this important result can be found in Sec.|4] 

Appendix B: Radio observables, synchrotron radiation and Stokes parameters 

This section introduces the notation to describe radio observables. For our statistical approach, synchrotron radiation is the funda- 
mental observed quantity on which our deduction is based. Since we are attempting to infer properties of the magnetic field statistics 
through statistics of the radio synchrotron observables, we need a clear and compact notation of these observables. 

All accelerated charges emit electromagnetic radiation. If accelerated by a magnetic field, the radiation is called cyclotron 
radiation in case of nonrelativistic and synchrotron radiation in case of relativistic velocities. With regard to astrophysics the latter 
is far more important. This is because of the much higher powe r radiated by relativistic p articles, since the total emitted power of an 
accelereated charge depends on y 2 , its Lorentz factor squared ( |Rybicki & Lightman|l979| l: 

P = lo-rCjSY^, (B.l) 

where f3 — v/c. Synchrotron radiation has a characteristic polarisation, with a high percentage of linear polarisation. Furthermore, 
the relativistic beaming effect confines the energy radiated within a cone around the direction of the moving charge. Due to these 
features and since a large fraction of astrophysical synchrotron emission falls into radio wavebands, it is relatively easy to detect 
and provides us with an excellent way to observe and study magnetic fields. 

Synchrotron radiation is mainly emitted by relativistic electrons. Other charged particles like protons contribute far less to the 
radiated power due to their larger mass. Following Rybicki & Ligh tman| ( |1979[ ) we assume a power-law distribution of the cosmic 
ray electron energies with spectral index p 

N(y) dy = Cy- p dy, (B.2) 

where C is a normalisation factor which determines the number density of relativistic electrons. The total power emitted per unit 
volume and unit frequency by such distributed electrons, assumed to have an isotropic pitch-angle distribution, is then given by the 
integral over N(y) dy times the single particle radiation spectrum. It can be shown that this leads to a power law in synchrotron 
emmissivity (Rybi ckT& Lightman|1979| : 

jccaj- — B ± 2 C. (B.3) 

In the end we are interested in observable quantities, namely the total and polarised intensity of the observed region. In terms of 
emissivity j they are: 



total intensity /(xj_) = I dz j(x), and (B.4) 
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polansedintenstty P(x x) = j dz jW( p) exp(2*(x)), (B.5) 

where f(p) = (p + l)/(p+7 /3) is the polarisation fraction and the integrals are along the line of sight from the source to the observer. 
The angle x is the polarisation angle of the radiation. With x we can introduce the effect of Faraday-rotation into our formulas. 
Faraday rotation is the rotation of the polarisation plane of a linearly polarised wave in a medium with a non-scalar dielectric 
constant due to a magnetic field. In such an environment, the dielectric constant differs for left and right circular polarisation 
( [Rybicki & Lightman|1979) : 



0)1 

g=l- , ' C . (B.6) 
a>{a> + u>b) 



where to denotes the frequency of the wave and a) p - eB/mc is the cyclotron frequency. Plus and minus signs denote the case for 
right and left circular polarisation respectively. If we describe a linearly polarised wave as a superposition of a wave with right 
circular polarisation and a wave with left circular polarisation, the linear polarisation plane will not remain constant. 
Including Faraday rotation, the total angle x at the location of the observer is given by 



Ar(x)=Aro(x) + A 2 ^(x). (B.7) 



Hereto denotes the polarisation angle at the position of emission whereas the Faraday depth 0(x) is defined as (Kronberg et al. 
2008) 



e 3 r 

^ ~ t TJ dzn e B 3 . (B.8) 

Z*7T 11 e C J 



It describes the phase angle through which the electric vector rotates due to Faraday rotation. For this study, it is assumed that 
the electron density n e is constant in order to simplify the c alculations. Furthermore, we restrict ourselves to observations at short 



wavelength A, thus ignoring the term exp(2iA 2 <f>(x)) in ( |B.5[ ), which would lead to so called intrinsic Faraday-rotation. Nevertheless 
4> is used as an independent radio observable, allowing us to infer information about the component of the magnetic field parallel 
to the line of sight. In real observations this would be possible by finding a background source to probe <p in the observed region 
(see Fig. [TJ. Sometimes we refer to the notion rotation measure (RM) instead of Faraday depth. Correctly this describes the factor 
between the rotation angle and A 2 obtained through observations. However, we use it in the context of real observations, because it 
is the more conventional term than Faraday depth when dealing with data. 



7(x ± ) and P(x ± ) can be expressed most suitable using the Stokes parameters. After some calculations (see e.g. |Rybicki & 
Light man|1979| Waelkens et al. 2009 ) the first three Stokes parameters I,Q and U can be determined for synchrotron radiation of a 



power-law spectrum distributed, isotropic, relativistic electron population : 

/ = 2F(p)w S r j dziBj + Bl^iBl + Bj), (B.9) 
6 = 2G(/))w a ? ! j dz(B\+Blf^(B\-B\), (B.10) 

U = 2G(p)o)^ j dz (B? + S|)^2BiB 2 - (B.ll) 

The two functions F(p) and G(p) are expressed in terms of physical constants and gamma functions of the spectral index p (see 
IWaelkens et al.|2009||KiBiin & Biermann|1998| l: 

W-m^CIXj-iwf*^. (B.12, 
5LTi L m e c L 3e 4 12 4 12 p + 1 

G( P ) = -^(^^Cixf - -Uxf + (B.13) 
32n £ m e c £ 3e 4 12 4 12 

Throughout this study, the value p — 3 has been adopted for the spectral index of the electron distribution. This not only simplifies 
calculations, it is also a reasonable choice from a physical point of view. Typical values for the spectral index of relativistic electrons 



in our galaxy measured directly by cosmic rays on Earth or indirectly via their synchrotron emission are around p a 2.7 ( Amsler 
|et al.|2008] >. Deviations from p — 3 can be added to the results of this work later in terms of corrections. 

Finally we can state /, P and <p in the form to be used in this study. During our calculations, all fore factors will be suppressed 
for convenience, resulting in 

1 = J dz(B\ + B\), (B.14) 
P = j dz{B\-Bl + 2iBiB 2 ), and (B.15) 
= j dzB^ (B.16) 
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Appendix C: Identity of the covariance matrix with the correlation tensor for Gaussian statistics 



In Sec. |2.3| we use the fact, that the covariance matrix of a Gaussian probability distribution of the magnetic field is identical to the 
magnetic correlation tensor. From the Gaussian probability distribution 



the identity is easily to shown: 



g(B,M) = exp[--B + M l B] 

V|2ttM| 2 



M = (BB 1 ) = j DB @{B, M)BB^ 

-(SSF)L-*5**1 



= ([M + MJ] exp[^y + M7)]| y=0 = M 



In this case, we have used the generating-function technique which is explained in 2.3 



